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Abstract 



A major topic of (classical) ergodic theory is to examine qualitatively how the 
phase space of dynamical systems is penetrated by the orbits of their dynamics. We 
CD ■ consider interacting qubit systems with dynamics according to 4-local Hamiltonians 

and continuous quantum random walks. For these systems one could use the von 
Neumann entropy of the time-average to characterize the mixing properties of the 
pi • corresponding orbits, i.e., what portion of the state space and how uniformly it is filled 

out by the orbits. We show that the problem of estimating this entropy is PSPACE- 
C [ hard. 

1 Introduction 

^ ' Ergodic theory studies the long-term average behavior of dynamical systems. The 

■ continuous time-evolution is described by a one-parameter family {Tt : t G M} of maps 

of the phase space X into itself. A major topic of classical ergodic theory is to examine 
qualitatively how the orbits {Tfx} (where x is the initial state) penetrate the phase 
space X. If the orbits "spread uniformly" over X, then the system is called ergodic. 
Ergodicity of dynamical systems plays an important role in statistical mechanics and 
thermodynamics because in the ergodic case the time averages and ensemble averages 
of physical variables are equal (see jl] for a rigorous study of these questions) . 

In this paper we consider the case that X is the Hilbert space of many qubits, 
and {Tt : t G M} is the unitary one-parameter group {e~*^* : t G M} generated by 
the system Hamiltonian H. We study the computational complexity of estimating the 
von Neumann entropy generated by Hamiltonian time evolutions if we average over 
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the time. The motivation for considering the entropy of the time average is that it 
quantifies in a certain sense what portion of the state space and how uniformly it 
is fihed out by the trajectory according to the time evolution. A second interesting 
property of the entropy of the time-average is that it could be used to quantify how 
far a pure state is from all eigenvectors of a Hamiltonian (this question was raised 
in The reason is that the entropy of the time-average of the initial pure state 

|$)($| can only be large if the initial state is a superposition of many eigenvectors 
(corresponding to different eigenvalues) of H with equal amplitudes. We explain this 
later in more detail. 

The time-average is defined as follows: 

Definition 1 (Time average) 

Let H be the Hamiltonian of a quantum system TC. For an arbitrary state p on TC its 
time average is the state p defined by 



It is well-known that the time average depends on the spectral decomposition of the 
Hamiltonian (a proof of the following lemma is included in the appendix for complete- 
ness). 

Lemma 1 The time-average p of p is given by 



where Px are the orthogonal projections onto the eigenspaces of H for different eigen- 
values A. 

Before we can formulate question on the complexity of estimating certain properties of 
time averages we have to restrict the class of considered Hamiltonians. This is because 
physically realistic Hamiltonians must satisfy some locality condition. One way to 
formalize this is the following definition. 

Definition 2 (local operator) 

Let TC := TCi • • • <^TCn be the tensor product Hilbert space of n Hilbert spaces. We call 
an operator k-local if it is a sum of operators that act on at most k tensor components 
non-trivially. 

If each tensor component is the state space of a physical particle, then the natural 
interactions are usually pair- interactions, that is, 2-local operators in the sense of the 
above definition. 

Nevertheless, for the Hilbert space Ti := (C^)*^" of n interacting qubits it may also 
be physically reasonable to allow A;-local interaction for k > 2. One reason is that they 
may describe effective Hamiltonians. A second reason is that there is not necessarily a 
one-to-one correspondence between qubits and physical particles (for instance, several 
qubits could describe the state of one particle) . 




(1) 



p = Y^p^pp^ 



A 
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We consider time-evolutions according to 4-local Hamiltonians and continuous quan- 
tum random walks. The time evolution e~*-^* is called a continuous quantum random 
walk if H is the adjacency matrix of some graph 8^. Recall that a matrix A is called 
an adjacency matrix if it is symmetric, has only Os and Is as entries, and all its entries 
on the diagonal are 0. Furthermore, we say that a matrix is a /c-local adjacency matrix 
if it is a sum of A:-local operators that are adjacency matrices themselves. We will 
consider 5-local adjacency matrices. 

We show that the problem estimating the von-Neumann entropy of time-averages of 
computational basis states, where the time-evolution is described by 4-local Hamiltoni- 
ans and 5-local quantum random walks, is PSPACE hard. The problem of estimating 
the density matrix of the time-average reduced to a qubit is also PSPACE-hard. 

Due to Lemma 1 it is clear that the entropy generation depends on the whole 
spectrum of the Hamiltonians. Therefore, these results should be compared to the 
results on the complexity of estimating the minimal eigenvalue (ground state energy) of 
local Hamiltonians. The problem of determining the lowest energy value of a (classical) 
spin-spin interaction of Ising type is known to be NP-complete For interacting 

qubits the problem of determining the lowest energy value is even QMA-complete 
( "Quantum-NP" ) if one allows 3- local interactions only |H1 [7] . Recall that in these NP 
and Quantum-NP problems the task is not to determine the lowest eigenvalues with 
high precision. The demanded accuracy is only inverse polynomially in the number n 
of interacting qubits. 

Note that the entropy of the time-average of a pure state p := |^)($| could be used 
to measure how "far away" |$) is from all eigenstates of H. This question was studied 
in the authors constructed quantum states that are far away from any eigenstate 
of any non-trivial local Hamiltonian, in the sense that |||<I>) — l^*)!! is greater than some 
constant lower bound for all eigenstates |^') of all fc-local (for a fixed k) Hamiltonians, 
independent of the form of the Hamiltonians. The measure for the distance of a pure 
quantum state |$) from the spectrum of a Hamiltonian used in |2] is 

d{\(b),H) := min min |||<I>) - |^a)|| , 
{I-i-a}} a 

where {|\I'a)} runs over all sets consisting of eigenvectors of H corresponding to different 
eigenvalues A. One easily verifies that the distance is equivalent to 

min min2(l — IoaP) , 
{|*A>} A ' 

where \^) = ^\0t\\^ \) . Therefore, the bound d{\^),H) is maximal if and only if 
\^) is an equally weighted superposition of eigenvectors corresponding to different 
eigenvalues of H. Using Lemma 1 one see that this is equivalent that the entropy of 
the time average of 1$) is maximal; the maximally possible entropy is the logarithm of 
the number of different eigenvalues of H. Based on these observations, our results may 
be also interpreted as a proof of the computational difficulty of deciding for a given 
quantum state how far it is from all eigenstates of a local Hamiltonian. 
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2 Stating the problem 



The techniques we use here to prove PSPACE-hardness of estimating the entropy of 
time-averages are closed related to those in [3] where we studied the complexity of 
measuring local observables. The general idea is to describe the computational steps 
carried out by a (reversible) space-bounded Turing machine by a quantum circuit 
consisting of elementary quantum gates or Toffoli gates and then to encode the circuit 
into a local Hamiltonian (this will be explained in Sections 3 and 4). 

Theorem 1 (PSPACE-hardness) 

Let (Ai) be a uniformly generated family of A-local Hamiltonians (or 5 -local adjacency 
matrices) on the Hilbert spaces Hi consisting of poly (l) qubits. Then the problems of 
(1) estimating the von-Neumann entropy of the time-average S{p) and (2) estimating 
the reduced state of p to single qubits are PSPACE-hard. 

More precisely: (1) Deciding if the entropy of the time-average of a computational 
basis state ofTCi according the dynamics given by Ai is smaller than ai or greater than 
bi, where bi — ai > poly{l), is PSPACE-hard. 

Let pi denote the reduced state of the time-average p to the ith qubit ofTCi. (2) 
Deciding if {l\pi\l) < 1/di or {l\pi\l) > 1/2 - 1/di, where di = poly {I), is PSPACE- 
hard. 

We prove this theorem in Section 4. The proof is based on a characterization of the 
complexity class PSPACE with quantum circuits that is explained in the next section. 



3 Characterizing PSPACE by circuits 

The complexity class PSPACE is usually defined with respect to the Turing machine 
model {Oj. PSPACE is the class of all languages recognizable by polynomial space 
bounded deterministic Turing machines that halt on all inputs 10 . 

In the following we work with a characterization of PSPACE with respect to circuits 
consisting of elementary quantum gates (Toffoli gates only). This characterization is 
proved in Theorem 1 in 3 . 

Theorem 2 (PSPACE) 

For every language L in PSPACE there is a polynomial-time uniformly generated family 
of quantum circuits (V;)/gN. Here I denotes the length of the input x. Each Vi consists 
of si = poly{l) elementary quantum gates (or Toffoli gates^) and acts on mi = poly{l) 
many qubits. The circuit Vi decides whether an input string x is an element of L in 
the following sense. 

There is a polynomial-time computable natural number ri such that the ri-fold con- 
catenation ofVi solves the corresponding PSPACE problem, i.e. 

V[Wx) ® |y) ® |00 . . . 0)) = \x) ®\y® f{x)) |00 . . . 0) , 

^In we have used elementary quantum gates instead of Toffoli gates. But it is clear that the transitions 
performed reversible Turing machines can be simulated by a circuit consisting of Toffoli gates only. This is 
because the Toffoli gate is universal for reversible computation. The reason why we use the Toffoli gates is 
that it is possible to construct from such a circuit a Hamiltonian is an adjacency matrix. 
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Figure 1: Schematic representation of how V works 



where f is the characteristic function of L. That is f{x) = 1 if x G L and f{x) = 
otherwise. The vector \x) is the basis state given by the binary word x G {0,1}', the 
vector \y) is the state of the output qubit and |00 . . . 0) is the initial state of mi — I — 1 
ancilla qubits. 

The intuition behind this construction is as follows. Recall that irreversible TMs (used 
to define PSPACE) can be simulated space-efficiently by reversible TMs |1H I12j. The 
computational steps of a reversible TM can be simulated efficiently by a quantum circuit 
U. This circuit acts on registers representing the state of the head of TM, the current 
tape position, and a polynomial portion of the tape (see jSJ for an explicit construction 
of the quantum circuit U) and realizes permutations of basis states corresponding to 
the transitions of the TM. By augmenting U with some control logic, we can achieve 
that once the TM moves into a final state, that the new circuit performs some idle 
cycles, flips the output qubit if and only if the answer is yes after exactly r//2 steps, 
then performs some idle cycles, and reverses the computation. This is schematically 
illustrated in Figure ^ It is necessary to perform idle cycles because we sometime 
do not know exactly after how many steps the TM moves into a final step, but only 
know a lower bound. By including the idle steps we guarantee the computational time 
depends only on the input length. 

Remark 1 We can consider an output register instead of an output qubit. It is also 
possible to construct a quantum circuit such that 

Vp{\x) (g) \y) (g) |00. ..0)) = \x) \y + f{x) mod 2") O |00. ..0) , 

for any c that is polynomial in I. 

4 Constructing the Hamiltonians / 
the quantum random walks 

Starting from the family (V/) of circuits we construct a uniformly generated family of 
4-local Hamiltonians (5-local adjacency matrices) (^i)igN such that the time-average 
of the state 

|x) |0) ® I00...0) I100...0) 



5 



and its entropy encodes the solution the instance x of a PSPACE problem; the first 
three tensor components are as in Theorem 2 and the last one corresponds to the 
new register clock that is needed in our construction. The construction is based 
on Feynman's construction of a computer whose dynamics is an autonomous time- 
evolution [13 EI. 

To explain the construction of the Hamiltonian (adjacency matrix A) let y be a 
circuit as in Theorem 2 and s be its size s, that is, the number of elementary gates 
(Toffoli gates). We need a register clock indicating which gate is applied. It consists 
of s qubits. The allowed states of the register clock are of the form |0 • • • 010 • • • 0) 
indicating which gate of V is applied currently. We denote by the elementary 
quantum gates (Toffoli gates) of V. 

We first define the forward-time operator 

F = To® |l)2(0|i® |0)i(l|o + 

Tl® |1)3(0|2® |0)2(l|l + 

T,_i® |l)i(0|o® |0),(l|,_i, 

where the operators |0)fc(l|fc and |l)fc(0|fc are the annihilation a and creation operators 
, respectively, acting on the ith qubit of the clock. The backward-time operator is 
defined as the adjoint of F. The Hamiltonian is defined as 

A:= F + F'^ . 

Now we show that if V consists of Toffoli gates only, then A is a 5-local adjacency 
matrix of some graph in the computational basis, that is, A is symmetric, has only Os 
and Is as entries, and its entries on the diagonal are all Os. Since the matrix describing 
the action of the Toffoli gate (in the computational basis) contains only Os and Is, the 
sums 

Tk ® |l)fc(0|fc ® |0)fc+i(lU.+i + Tk \0)k{l\k |l)fc+i(OU+i 

are adjacency matrices that are 5-local operators. A is the sum of the above matrices. 
Therefore, it is symmetric and has only Os on the diagonal. It remains to show that 
no Is of these matrices can meet when summing. This is done by checking that the 
operators 

O-k^k+l mod s , O-kO-k+l mod s ' 

lead to orthogonal states for all k = 0, . . . , s — 1 when applied to any computational 
basis state. We see this by observing that \i k — I mod s > 2 then we have 

{h\a\ak+i aia\^^\h) = and {b\aka\^-^ aia\^^\h) = 

because the operators a\ak+i aia\j^-^ and aka\j^-^aia\j^^ either map |6) onto the zero 
vector or change the bits of |6) at the non-overlapping positions {k, k + 1) and {I, ^ -|- 1). 
Similarly, if / = A; -|- 1 then we have 

{h\a\ak+i ak+ia\j^^\h) = and {b\aka\^^ ak+ia\^^\h) = 
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because aj[.afc_|_i ak+ia]^_^_2 is the zero operator and 0^0^.^^ afe+iaj[.^2 either maps \b) onto 
the zero vector or changes the bits of |6) at the non-overlapping positions k and k + 2. 

Therefore, only computational basis states appear in A\b) with coefficients all equal 
to 1. Consequently, we have either = or = 1. This show that A is a 

5-local adjacency matrix if V contains Toffoli gates only. 

We denote the linear span of the vectors 

F^\-^o) forjGN 

where 1^*0) := \xOQ - ■ - 0) |100 • • • 0) as C Here the first part of the tensor product 
denotes the register where V acts on (input and ancilla registers) and the second 
component is the register clock. All states of this orbit are orthogonal until one has a 
recurrence to the initial state |^'o)- This can be seen as follows: If the register clock 
is in an allowed state there is only one summand of F that is relevant. Its action on 
the clock is simple since it moves the symbol 1 to the next qubit (this dynamics of the 
clock may be interpreted as a propagation of a spin- wave). Therefore it is clear that 
the first s — 1 states are orthogonal. The whole circuit y is a classical logical operation 
which permutes basis states. Therefore the state F*|^'o) is either orthogonal to |5'o) 
or both states coincide. Along the same line we argue that all states of the orbit are 
orthogonal until a state coincides with the initial state. Hence F acts as a cyclic shift 
on O. 

For our construction is is essential that the dimension of the orbit depends on the 
solution of the PSPACE problem. It is 2sr if f{x) = 1 and sr if f{x) = 0, where f{x) 
is the answer for the instance x of the considered PSPACE problem. We denote the 
dimension of the orbit O by d. Let be a primitive complex d-th root of unity. The 
eigenvalues of F restricted to O are 

LO ,UJ , 0^2, . . . ,iO 

Furthermore, |^'o) is a superposition of all eigenvectors of F restricted to O with equal 
amplitudes. All this follows from properties of the cyclic shift operator. Since F and 
commute on O the eigenvalues of A restricted to O are oj^ + oj^ = 2 cos(27rj/c/) for 
/c = 0, ...,d-l. 

Set p := |^'o)(^o| and := 2 cos{2TTk / d) . Note that d is always even; this is 
because r in Theorem 2 is even. Therefore, there are two eigenvalues with multiplicity 

1 (corresponding to the real eigenvalues +1 and — 1 of F) and d — 2 eigenvalues with 
multiplicity 2 (corresponding to the complex eigenvalues of F). We have 

fc,;=o 

where \k) denotes the kih. eigenvalue of F restricted to O. For the time-average we 
obtain with Lemma 1 

d-1 



p = \Y.'^>^k = k\\k){i\^ 



k,l=0 
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where [A^ = A;] is 1 if both eigenvalues are equal and otherwise. If we permute 
the order of the eigenvectors \k) such that eigenvectors corresponding to the same 
eigenvalues are adjacent, then p looks in this basis as follows: 



/ 1 



1 1 
1 1 



v 



1 1 

1 1 / 



(2) 



Using the form of the time-average in , its entropy is easily computed 

S{p) = log2 d-id- 2)/d « log2 d-1. 

Therefore, the difference in entropy of the time-average of the computational basis 
state l^o) encoding a "yes" -instance and a "no"-instance is at least c/2 bits, where c 
is the size of output register (see Remark 1). This proves the PSPACE-hardness of 
estimating the entropy of the time-average as stated in Theorem 2. 

Let the output register consist of only one qubit. Obviously, the reduced state of 
time-average to the output qubit is |0)(0| for a "no" -instance. For a "yes" -instance we 
know that the output qubit of the states := F-' l^'o) is 1 for j = d/4, . . . , 3d/4 — 1 
and for all other j. To see this, recall how V works: at step r//2 the circuit V changes 
the output qubit from to 1, at step r/ we obtain the initial state with the output 
flipped, at step r; -|- ri/2 the output is flipped from 1 to 0, and finally in step 2r; we 
obtain the original input state. 

Therefore, the occupation probability of the level |1) of po^t {p restricted to the 
output qubit) is exactly 

3d/4-l 
j=d/4 

Let us denote by -P(j) the probability to observe the state \^ j) if we measure the 
time-average in the computational basis. The probability P{j) is given by 

p{j) = {^^m,) 



1 1 



This seen by observing that \^ j) = Z]fc=o '^"'^l^)- Now we make use of the Diaconis- 
Shahshahani bound (see Appendix) applied to the cyclic group Z^. The characters of 
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1^^ are given by Xa(j) = ^""^ ■ The coefHcients of the Fourier transform of the probabihty 
density P = (P(0), . . . , P(d - 1)) are 



3=0 



fe^/ ■.\k=\i j=o 

Note that the possible values for k — I are only the even numbers 2,4, . . . ,d/2 since 
we sum over all pairs {k,l) satisfying k ^ I and the condition cos(27rfe/(i) = cos{2nl/d) 
being equivalent to k = —j. 

If a is odd, then k — l + a^O mod d. But if a is even, then there is exactly one pair 
of {k, I) satisfying above properties and k — I + a = mod d. With these observations 
we obtain 

. Jo if a odd , . 

P{a) = < I .f (3) 
^ ^ [2 if a even 

Now we obtain with Diaconis-Shahshani bound a lower bound on the total variation 
distance between P and the uniform distribution U on the orbit states F-' |^'o) 

||F-l/|||v<iw2-l)J,<i^ 

Taking the square root and multiplying both by 2 we obtain 

d-i 

Using the triangle inequality the above bound implies that 

1111 

(l|Pout|l) = (l|p|l) = Yl ^ 



,,,, - v2ci'2 ^/2d 

J=d/4 

(Note that this is also true if we sum over any set containing half of all possible j.) 

In summary, we have proved that for a "yes" -instance the occupation probability 
(l|pout|l) is almost 1/2, whereas it is for a "no" -instance. 



5 Appendix 

Proof. Lemma 1: Let {|A;)} be a basis consisting of eigenvectors of H with corre- 
sponding eigenvalues A^. Express the state p in this basis 

P = J2^ki\k){l\. 

k,l 
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For the time-average we obtain 

-p = lim 1 r ^e-'('''-'^^'jki\k){l\ 



k,l=0 
d-l 



k,l=0 



where [A^ = A;] is 1 if both eigenvalues are equal and otherwise, A are the differ- 
ent eigenvalues of H, and Px are the orthogonal projection onto the corresponding 
eigenspaces. □ 



The proof of the lemma can be found in JS] . 
Lemma 2 (Diaconis-Shahshahani bound) 

Let A be an arbitrary abelian group, P an arbitrary probability distribution on A and U 
the uniform distribution on A. The Diaconis-Shahshahani bound on the total variation 
distance 

WP-Uhv ■■=IJ2\P(9)-Uig)\ 

g<^A 

is 

x¥=xo 

where the sum runs over all non-trivial irreducible characters of A and P denotes the 
Fourier transform of P, that is, P{x) = XlgeA x(o)-P(fi')- 
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